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ABSTRACT 

A critical discussion is presented of the data analysis applied by Bernstein, Freedman, & Madore 
(2002a, b) in their measurement of the Extragalactic Background Light. There are questionable as- 
sumptions in the analysis of the ground-based observations of the Zodiacal Light. The modeling of the 
Diffuse Galactic Light is based on an underestimated value of the dust column density along the line of 
sight. Comparison with the previously presented results from the same observations reveals a puzzling 
situation: in spite of a large difference in the atmospheric scattered light corrections the derived Extra- 
galactic Background Light values are exactly the same. The claim of the paper of a "detection of the 
Extragalactic Background Light" appears premature. 

Subject headings: diffuse radiation - techniques: photometric - dust, extinction 
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1. INTRODUCTION 

Bernstein, Freedman, & Madore (2002a,b), hereafter 
BFM02a, BFM02b, or BFM02 combined, have an- 
nounced the first detection of the Extragalactic Back- 
ground Light (EBL) from absolute photometry, with 
the mean values of 4.0(±2.5), 2.7(±1.4), and 2.2(±1.0) 
xlO^ergs^cm-^sr- 1 !- 1 at 3000, 5500, and 8000 A, 
respectively. The errors quoted are la uncertainties. Their 
method is based on the formula: 

^EBL = hot — Izh — IDGL, (1) 

where I t ot is the total sky surface brightness outside the 
atmosphere, Jzl is the Zodiacal Light (ZL), and Idgl is 
the Diffuse Galactic Light (DGL) surface brightness, all 
to be determined in the direction of the BFM02 target 
field at I = 206.6, b = —59.8 deg. Each one of the three 
components is derived in BFM02 with a different method, 
/tot is measured above the atmosphere with the Hubble 
Space Telescope (HST) using broad-band CCD photome- 
try; Izl is measured from the ground, with the 2.5-m du 
Pont telescope at the Las Campanas Observatory (LCO), 
using spectrophotometry; and /dgl is estimated by us- 
ing a model for the scattering of starlight by interstellar 
dust. A very demanding task for the BFM02 method is 
set by the requirement that, for each of the two telescopes 
with different properties and different observing methods, 
the measured flux, Izl or 7 t ot, has to be separately cal- 
ibrated to the same scale. Iebl is only a small fraction, 
a few per cent at most, of I to t and 2zl- Therefore, the 
BFM02 method crucially depends on whether or not the 
very high absolute accuracy, of <1 %, needed in the mea- 
surement, calibration, and scattered light corrections for 
/zl and separately for itot is achieved. 

In this paper a critical discussion will be presented of 
the calibration and the atmospheric corrections applied in 
BFM02b to the ground based measurement of the Zodiacal 
Light. The assumptions for estimating Idgl are critically 
reviewed. In addition, I point out a puzzling situation 
which emergerges from comparison between the results of 
BFM02 and of the widely cited previous presentations of 
the same observations. 



2. GROUND BASED MEASUREMENT OF THE ZODIACAL 

LIGHT 

Two atmospheric corrections have to be applied: (1) 
extinction, and (2) tropospheric scattered light. The ob- 
served night sky brightness, I ODS (A, i, A), towards the tar- 
get field is given by 



^obs (A, t, X) — 

/z L (A)e- T(A)x + [Jebl + J DG l](AK 

^agl(A,i, X) 



isca(A, t, X) 



r(\)X 



(2) 



where A is the wavelength, t the time of the observation, X 
the airmass, r(A) the atmospheric extinction coefficient for 
unit airmass, 7 sca the tropospheric scattered light, and 7 ag i 
the airglow as observed from the ground (including atmo- 
spheric attenuation and scattered airglow). The BFM02b 
method of separating the ZL from the airglow component 
is based on the assumptions that (1) the depths of the 
Fraunhofer lines in the spectrum of the Zodiacal Light are 
identical to those in the solar spectrum, and (2) the air- 
glow spectrum is uncorrelated with the solar Fraunhofer 
spectrum. 

The tropospheric scattered light, I sca (X,t, X), is the 
main obstacle in conducting accurate absolute diffuse sky 
photometry from the ground. Unlike the photometry of 
stars or small extended sources, no differential ON/OFF 
measurements are possible, and one must calculate the 
scattered light contribution coming from all the light 
sources above the horizon. The scattered light components 
which contribute to the "ZL-like" (Fraunhofer spectrum) 
signal are due to the all-sky distributions of the ZL itself, 
the Integrated Starlight (ISL), and the DGL. Each one of 
these components has the Rayleigh (R) and the Mie or 
aerosol (M) scattering part: 



J Bca (A,t,X) = l5 a (ZL) + /^(ZL) 
/* a (ISL) + /™ a (ISL) + JiL(DGL) - 



^a(DGL) 
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3. COMPARISON BETWEEN BFM02 AND THE PREVIOUS 
RESULTS OF BERNSTEIN ET AL. 

The results of BFM02 have been previously presented 
in in Bernstein (1998, 1999a,b, 2001), and distributed in 
a preprint form for some time (Bernstein, Freedman, & 
Madore 1999, 2000); in the following they are collectively 
referred to as BFM98-01. These results have been widely 
cited in the literature as the "EBL standard reference val- 
ues", see e.g. Barger et al. (2001); Hauser (2001); Jimenez 
& Kashlinsky (1999); Longair (2001); Madau & Pozzetti 

(2000) ; Pagcl (2002); Peebles (2001); Pozzetti & Madau 

(2001) ; Primack et al. (1999); Renault et al. (2001); Wright 
(2001). It is therefore important to point out the follow- 
ing puzzling situation which results from the comparison 
between the BFM98-01 and BFM02 results: 

In BFM98-01, the reduction method for the Zodiacal 
Light measurements differed in a fundamental way from 
that in BFM02b. Instead of Eq.(2), the following formula 
was applied: 

4>bs(A,t,A) = 

hL(X)e- TWX + [Jebl + /dglPK^ + 
4gi(A,*,X) (3) 

i.e. the atmospheric scattered light term 7 sca (A,t, X) was 
completely omitted without explanation, thus neglecting 
this fundamental aspect of the diffuse night sky photom- 
etry. Only in BFM02b has an Appendix on scattered 
light model calculations now been added. Using the re- 
sults of these calculations a lower limit to I sca at 4650 
A can be estimated by adopting X = 1.1 for the air- 
mass, i.e. a value at the lower end of the BFM02b airmass 
range. 7 sca amounts to — 18% of Izl- It consists of the 
Rayleigh (—10% of Izl) and Mie (—4% of Izl) components 
of scattered Zodiacal Light, as well as of a scattered ISL 
component with a "ZL-like" spectrum (—4% of Izl) (see 
Sect. 4 and 5, and Figs. 20, 21, and 30 of BFM02b). This 
value translates, for the airmass of 1.1, to an outside-the- 
atmosphere value of —22% of Izl which, for the BFM02b 
value of Izl = 109.4xl0 _9 crgs _1 cm _2 sr _1 A _1 at 4650 A , 
corresponds to ^24x 10 _9 ergs _1 cm _2 sr _1 A _1 . 

Since not subtracted, I sca was erroneously included into 
the Zodiacal Light, 7 Z l, in the BFM98-01 analysis. Thus, 
the difference in the reduction methods necessarily should 
have resulted in a —22% larger ZL value in BFM98-01 than 
in BFM02b. This -22% difference in I Z l is at least 7 times 
as large as the EBL value of BFM02 at 5500 and 8000 
A. However, in BFM98-01 the EBL intensities at 3000, 
5500, and 8000 A are identical to the values in BFM02 
of 4.0, 2.7, and 2.2x lO^ergs^cm^sr^A- 1 to within 
^ 0.1 x 10 _9 crgs _1 cm _2 sr _1 A _1 . It is unclear how this 
puzzling situation should be understood. 

4. SCATTERED LIGHT CORRECTIONS 

In an Appendix of BFM02b model calculations for the 
atmospheric scattered light are now presented. How- 
ever, the nature of the problem dictates that it is 
hardly possible to achieve an absolute accuracy of 
<lxl0 _9 crgs _1 cm -2 sr _1 A _1 as required in the BFM02 
method. Major problems are caused by e.g. the vary- 
ing properties of the atmospheric aerosols and ground re- 
flectance, as well as by the insufficiently known intensity 



distributions of the main light sources, the ZL, ISL, and 
DGL over the sky. 

The Moon, if above or up to a few degrees below the 
horizon, will give rise to a substantial atmospheric scat- 
tered light component (Krisciunas & Schacfer 1991). The 
observations at LCO, if carried out on 1995 November 27 
and 29 as announced by BFM98-01 and BFM02b, would 
have included time slots when the Moon was above the 
horizon. In an erratum (Bernstein, Freedman, & Madore 
2003) the authors have now removed this problem by stat- 
ing that their published dates were incorrect, the correct 
dates of their observing nights being 1995 November 24/25 
and November 26/27. 

4.1. Aerosol scattering 

BFM02b have adopted for the albedo of the aerosol par- 
ticles a value of a = 0.59. This value, given in Staude 
(1975), was calculated for an ad hoc particle composition 
with a refractive index of m = 1.5 — O.li. For a realis- 
tic aerosol composition, according to Garstang (1991) and 
McClatchey et al. (1978), the albedo is -0.94. This is 
a representative value over the whole wavelength range, 
3900 - 5100 A , of the BFM02b Zodiacal Light measure- 
ment and for a variety of different measured aerosol pop- 
ulations and conditions. The aerosol scattering compo- 
nent, I™ a = I™ a (ZL) + J s ^ a (ISL), has to be corrected for 
the effect of increased albedo by multiplying it with the 
albedo ratio 0.94/0.59. In order to obtain the contribu- 
tion to the extraterrestrial Izl value one has to multiply 
^sca by the extinction correction factor e T ^ x . 

According to BFM02b Appendix, I s ^ a (ZL) is 3.6 - 6.1 % 
of Izl at their minimum airmass of 1.02 and 5.2 - 8.1 % 
of Izl at their maximum airmass of 1.82. The lower value 
stands for A = 4920 A and the higher one for A = 3960 
A. ( Four selected wavelength slots at -3960, 4270, 4340, 
and 4920 A were used by BFM02b to derive their final Izl 
value at 4650 A.) 

41 (ISL) is not explicitly given in BFM02b Appendix. 
It is only a small fraction of the total (Rayleigh + Mie) 
scattered ISL which they estimate to be — 12 — 24 
xlO _9 ergs _1 cm _2 sr _1 A _1 . A rough estimate is obtained 
from 

C(ISL) « o(l - e-™*)/isL, (4) 

where tm is the aerosol extinction per airmass and 4sl ~ 
40 xlO _9 crgs _1 cm _2 sr _1 A _1 is the ISL intensity in the 
direction of the target. For the wavelength and airmass 
range, and the albedo a = 0.59, as adopted by BFM02b, 
4^ (ISL) « (0.9-2.4) xlO- 9 ergs- 1 cm- 2 sr- 1 A- 1 . A 
third of 41 (ISL) is estimated to have a solar type Fraun- 
hofcr spectrum and thus adds —0.3 - 0.8 % of Izl to the 
scattered ZL value. 

The average of I^ a e T ^ x for the 16 different aimasses 
and four wavelength slots of the BFM02b measure- 
ment gives an outside-the-atmosphere value of —7.6 % 
of Izl- With the aerosol albedo of 0.94 instead of 
0.59, this value should be scaled up to —12.1 % of 
Izl- The difference of —4.5 % of Izl corresponds to 
4.9xlO~ 9 ergs~ 1 cm~ 2 sr~ 1 A~ 1 at 4650 A which, when not 
subtracted, will artificially increase the Izl value (sec Ta- 
ble 1). 
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4.2. Ground reflection 

BFM02b have neglected the reflection from the ground. 
Light reflected from the ground is scattered a second time 
by the molecules in the atmosphere, back into the ob- 
server's line of sight. Because of the strongly forward scat- 
tering phase function of the aerosols, only the Raylcigh 
scattering by molecules is important here. The influence 
of the ground reflection on the Rayleigh scattered light 
intensity can be estimated using the tables of Ashburn 
(1954). The ground reflectance value obtained from the 
NASA MODIS/Tcrra Surface Reflectance database 1 is 
~8 %. This value is for a 100x100 km 2 area centered at 
the Las Campanas Observatory, for the same season of the 
year (November-December 2000 and 2001) as the BFM02b 
observations, and for a wavelength band (4590 - 4790 A , 
MODIS/Tcrra band 3) which closely matches the one used 
in BFM02b. According to the Ashburn (1954) tables the 
Raylcigh scattered light intensity, for the extinction range 
of BFM02b, has to be increased by 7.2% relative to the 
case of zero ground reflectance. 

According to BFM02b Appendix, J* a (ZL) is 8-18 % 
of Izl at their minimum airmass of 1.02 and 15-32 % of 
Jzl at their maximum airmass of 1.82. The lower value 
stands for A = 4920 A and the higher one for A = 3960 A. 
The value of 7* a (ISL) is not explicitely given in BFM02b, 
but it can be estimated using the total (Raylcigh + Mie) 
scattered ISL which BFM02b estimate to be - 12 - 24 
xlO- 9 ergs- 1 cm- 2 sr- 1 A- 1 . After subtracting J^JISL) 
according to Sect. 4.1 and adding J^: a (DGL) according to 
Sect. 4.3 one obtains the range i£ a (ISL) + I* a (DGL) = 
(14-27) xlO _9 crgs _1 cm _2 sr _1 A _1 . A third of it is esti- 
mated to have a solar type Fraunhofer spectrum and thus 
adds ^4.6 - 8.3 % of Izl to the scattered ZL value. 

The average of the ground reflection correction referred 
to outside-the-atmosphere, AI^: a e r ^ x , at the 16 different 
aimasses and four wavelength slots of the BFM02b mea- 
surement, gives a value of ^2.1 % of Izl- This corresponds 
to 2.3 xlO _9 ergs~ 1 cm _2 sr _1 A~ 1 at 4650 A which, when 
not subtracted, will artificially increase the Izl value (see 
Table 1). 

4.3. Diffuse Galactic Light as source 

BFM02b have constructed an approximate model for the 
intensity distribution and spectrum of the ISL. They find 
that their averaged intensities are within 10 % of the star 
count integrations of Roach & Megill (1961). However, 
the ISL is only a part of the Galactic radiation field which 
contributes to the tropospheric scattered light, with the 
other part being the DGL. The contribution of the DGL 
to the total Galactic radiation field can be estimated us- 
ing the tabulation of the ratios Idgl/Iisl- These data 
are based on the Pioneer 10 measurements at 4400 A by 
Toller (1981) and are reproduced in Table 39 of Leinert 
et al. (1998). The Idgl/Iisl ratio varies as a function 
of galactic latitude, from ~0. 25-0. 34 at low galactic lati- 
tudes to ^0.12 at high latitudes. There is no pronounced 
dependence of this ratio on the galactic longitude. To es- 
timate the total contribution of the DGL to the Galactic 
radiation field a weighted mean value of the ratio over the 
galactic latitude range was formed. The overall ratio is 

1 Available at http://modis-land.gsfc.nasa.gov/mod09/ 



0.24, i.e. the pure ISL sky brightness used by BFM02b 
as the Galactic illumination source has to be increased 
by 24%. This 24% correction should be applied to the 
Raylcigh scattering component only. Because of strong 
forward scattering, the Mie (aerosol) component couples 
to the DGL in the viewing direction at high galactic lati- 
tudes where Idgl/Iisl ~ 0.12. 

According to BFM02b, the scattered (Rayleigh+Mie) 
ISL is 12 - 24xlO~ 9 crgs~ 1 cm~ 2 sr~ 1 A~ 1 over their 
airmass range of 1.02 to 1.82. The Mie com- 
ponent has been estimated in Sect. 4.1 to be 
7^ a (ISL) « (0.9 - 2.4)xlO- 9 crgs- 1 cm- 2 sr- 1 A- 1 . 
Thus, the Rayleigh component is i s ^ a (ISL) « (11 — 
22)xlO -9 ergs~ 1 cm~ 2 sr~ 1 A~ 1 . The atmospheric scat- 
tered DGL component is 24% of it, i.e. 7* a (DGL) « 
(2.6 — 5.2)xlO _9 ergs _1 cm _2 sr _1 A _1 . However, as dis- 
cussed in BFM02b, it is not the total scattered flux which 
is important, but rather the strength of the Fraunhofer 
spectral features which are in common with the Sun. Ac- 
cording to BFM02b, the strengths of the spectral features 
they use in their analysis are approximately 1.5 to 3.8 
times weaker in the DGL than in the ZL spectrum. To 
obtain a reasonable lower-end estimate for the "ZL-like" 
component I use here a factor of 3.0, leading to a range of 
(0.9-1.8) xlO- 9 ergs- 1 cm- 2 sr- 1 A- 1 . 

The average of / s ^ a (DGL)c T ( A ) x at the 16 different 
aimasses and four wavelength slots of the BFM02b 
measurement gives an outside-the-atmosphere value of 
^1.6xlO _9 ergs _1 cm _2 sr _1 A _1 at 4650 A which, when 
not subtracted, will artificially increase the Izl value (see 
Table 1). 

5. CALIBRATION OF THE ZODIACAL LIGHT 
MEASUREMENTS: APERTURE CORRECTION 

As stated in BFM02a,b the calibration of extended uni- 
form surface brightness photometry differs in an essen- 
tial way from point source photometry. The following two 
standard issues are in common with point source and sur- 
face photometry: 

(1) Point source calibration (including the atmospheric ex- 
tinction effects); and (2) Calibration of the fiducial stan- 
dard star system. 

However, in the calibration of extended uniform surface 
brightness photometry one requires knowledge of two ad- 
ditional aspects: 

(3) The solid angle subtended by the spectrometer or pho- 
tometer aperture or CCD detector pixel; and (4) The aper- 
ture correction factor which accounts for the loss of flux of 
the standard star outside the spectrometer or photometer 
aperture. 

Because of scattering by micro-roughness and dust par- 
ticle contamination on the optics, and atmospheric small- 
angle scatting, the point spread function of a telescope 
extends beyond the aperture size of 1-20 arcsec normally 
used in point source photometry The aperture correc- 
tion factor, T(A), is the fraction of the flux from a point 
source that is contained within the aperture. The fraction 
1 — T(A) of the point source flux is lost outside the aper- 
ture. In the measurement of a uniform extended source 
the situation is different: the flux which is lost from the 
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solid angle defined by the focal plane aperture is com- 
pensated for by the flux which is scattered and diffracted 
into the aperture from the sky outside of the solid an- 
gle of the aperture. Therefore, the intensity in units 
of ergs _1 cm _2 sr _1 A _1 of an extended uniform source is 
given by 

m . « C(4 ,5, 

where C(A) is the signal in instrumental units (DN s -1 ), 
Q the solid angle of the aperture in steradians, and S(X) 
the sensitivity function in units of erg s _1 cm~ 2 A -1 / (DN 
s _1 ) determined from the standard star observations. 

The aperture correction is an additional factor specific 
to the absolute surface photometry of extended uniform 
sources. It does not appear in point source photometry, 
nor in the photometry of small extended sources where the 
sky background can be measured next to the object. The 
aperture correction factor is not "automatically" taken 
care of by simply using in the standard star observations 
the same measuring arrangement as was used in generating 
the standard star system. 

The accuracy of the EBL measurement of BFM02 cru- 
cially depends on how consistently the surface brightnesses 
for the HST and du Pont telescope can be calibrated to 
the same scale in spite of the different properties of the two 
telescopes. In practice, as also pointed out by BFM02, the 
majority of error in calibrating a uniform surface bright- 
ness source comes from the accuracy with which the large- 
angle PSF is determined. In fact, the two telescopes differ 
strongly with respect to the aperture correction factor. 

The surface (spectro)photometry of BFM02b at LCO 
was calibrated by standard stars which were observed 
through a slit of 10.8 arcsec width. In order to compen- 
sate for the light lost outside of the slit, BFM02b mea- 
sured the PSF up to a radius of ~60 arcsec. They found 
that the aperture correction factor for a uniform-surface- 
brightness, aperture-filling source is T = 0.963. However, 
in their analysis, BFM02b did not take into account that, 
in order to include 100% of encircled energy into the star 
image, the integration must be extended far beyond 60 
arcsec from the axis. 

A widely used compilation of data for the profile of a 
star image was presented by King (1971). Outside the 
central disk of ~10 arcsec radius, the image shows a more 
slowly declining halo or aureole which is well represented 
by an inverse-square law of intensity over a factor of 1000 
in angular distance. Between 1 and 100 arcmin the aureole 
contains about 5% and between 1 arcmin and 5 deg about 
6% of the star's flux. Later work by e.g. Capaccioli & de 
Vaucouleurs (1983) (McDonald 0.9-m and 2-m telescopes), 
Surma, Seifert, & Bender (1990) (Calar Alto 1.23-m), 
Racine (1996) (CTIO 4-m), Mackie (1992) (KPNO 0.6-m 
Burrel Schmidt), Middlemass, Clegg & Walsh (1989) (La 
Palma ING 2.5-m), Piccirillo (1973) (Goethe Link 16-inch 
and McDonald 2-m), and Uson, Boughn, & Kuhn (1991) 
(KPNO No. 1 0.9-m) has demonstrated that very similar 
functional dependences (i.e. inverse square law), but in 
some cases substantially higher aureole energy fractions, 
are obtained for other telescopes equipped with photo- 
electric photometers, CCD cameras or spectrographs. For 
example, Capaccioli & de Vaucouleurs (1983) find an inte- 
grated aureole contribution of ~10% between 30 arcec and 



1.5 deg. Furthermore, in Uson, Boughn, & Kuhn (1991), 
the aureole stray radiation level is ~1.5 times higher than 
in King (1971). The probable reason for the aureole is 
light scattering by the imperfections of the telescope op- 
tics, such as microripples and dust on the mirrors (Beckers 
1995; Roddier 1995). 

Although no stellar image aureole measurements were 
available in BFM02b or BFM98-01 for the du Pont tele- 
scope, there are good reasons to assume that its charac- 
teristics are similar to the telescopes mentioned above. 
Including the aureole contribution to the stellar PSF, 
the aperture correction factor of BFM02b should be de- 
creased by 5 - 10%, from T = 0.963 to 0.86 - 0.91. 
The resulting outsidc-thc-atmosphere ZL value should 
thus be reduced by 5 - 10%, corresponding to 5.5 - 
llxlO- 9 ergs- 1 cm- 2 sr- 1 A- 1 at 4650 A (see Table 1). 

For their HST measurement BFM02a used the WFPC2 
default calibration and adopted the standard aperture cor- 
rections as given in Holtzman et al. (1995). For their filters 
used, they applied an aperture correction factor of 90 % in 
moving from a 0'.'5 radius aperture to an "infinite" aper- 
ture (6'.'0 radius). Large angle scattering measurements by 
Krist & Burrows (1994) at 20"- 60"have shown that the 
scattered light level in HST/WFPC2 is much below the au- 
reole levels of the ground-based telescopes as listed above. 
Uson, Boughn, & Kuhn (1991) and Surma, Seifert, & Ben- 
der (1990) give at 60"a stray radiation value of 15.8 mag 
arcsec -2 (normalisation to a th mag star), correspond- 
ing to 4.8xl0~ 7 arcsec -2 (normalisation to stellar flux = 
1), while the HST/MVFFC2 value is 7xl0~ 8 arcsec" 2 , i.e. 
a factor of ~7 lower. This is as expected since the HST 
optics are known to have extremely small microroughncss 
and are practically dust-free. Therefore, it is unlikely that 
any additional large-angle aureole correction to the HST 
aperture correction factor is needed beyond the 90% cor- 
rection applied by BFM02a. In any case it would be much 
smaller than the correction for the ground-based du Pont 
telescope. Thus, when taking the difference I to t ~ Izl, the 
error in the aperture correction factor for Izl (du Pont 
telescope) is not compensated for by a similar error in 7 t ot 
(HST) . 

6. MODEL ESTIMATES OF THE DIFFUSE GALACTIC LIGHT 

The BFM02a estimation of the DGL intensity is based 
on the 100 fxm surface brightness, Iiooj a $ extracted from 
the IRAS Sky Survey Atlas (ISSA). This Atlas gives the 
very low value of 0.4 MJy sr -1 for their target position. 
The authors have apparently not paid attention to the fact 
that the ISSA surface brightnesses cannot be utilised for 
absolute surface photometry; see The Explanatory Sup- 
plement to the IRAS Sky Survey Atlas, Wheelock et al. 
(1994), p. 1-6. From their I 100 estimate BFM02a derived 
the values A(H) = 0.47 x 10 20 c m - 2 and A v = 0.028 
mag. However, the value of N(R) can be directly ex- 
tracted from the Hartmann & Burton (1997) Atlas, care- 
fully corrected for stray radiation effects. It gives for the 
BFM02 target A(H) = 1.78 x 10 20 cm" 2 corresponding 
to Ay = 0.106 mag, i.e. 3.8 times as large as the value 
adopted in BFM02a. Using the COBE/DIRBE surface 
photometry to fix the zero point for the IRAS 100 /im 
data, Schlegel, Finkbeiner, & Davis (1998) produced an 
all-sky Galactic 100 fim. emission and optical extinction 
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Table 1 

Corrections to be applied to the data analysis of BFM02: A7 is the correction above-the- atmosphere. The 

CORRECTIONS 1-4 ARE TO BE SUBTRACTED FROM THE Izl VALUES, THE DGL CORRECTION (ITEM 5) IS TO BE ADDED TO THE 
IdGL VALUE AS GIVEN IN BFM02. THE LAST THREE COLUMNS GIVE THE RESULTING CORRECTIONS TO BE APPLIED TO THE 

BFM02 7ebl values at 3000, 5500, and 8000 A. In lower part of the table the sum of the corrections 1-5, the 
BFM02 7ebl values at 3000, 5500, and 8000 A , and the Jebl estimates of this paper after applying the 

corrections are given. 



Correction BFM02 This paper 


A/o a Applies to 


A/ebl" 

@3000 A @5500 A @8000 A 


1. Aerosol albedo 0.59 0.94 

2. Ground reflectance neglected 8% 

3. DGL as source of 

atmospheric scattering neglected included 

4. ZL calibration, 

aperture correction factor 0.963 0.86-0.91 

5. DGL modeling, 

line-of-sight A v m 028 m 081±0 m 025 


-4.9 ZL @ 4650 A 
-2.3 ZL @ 4650 A 

-1.6 ZL @ 4650 A 

-5.5to-ll ZL @ 4650 A 

+ 1.5to+1.9 DGL <§> all A 


+ 1.5 +4.7 +3.2 
+0.7 +2.2 +1.5 

+0.5 +1.5 +1.1 

+ 1.7-3.4 +5.3-10.6 +3.6-7.3 

-1 -1 -1 


Sum of corrections 1-5 




+3.4-5.1 +12.7-18.0 +8.4-12.1 


Jebl (BFM02) 
Jebl (this paper) 




4.0 + 2.5 2.7+1.4 2.2+1.0 
7.4-9.1 15.4-20.7 10.6-14.3 



a in units of 10 9 ergs 'cm 2 sr 1 A 1 

^Corrections in units of 10 _9 ergs -1 cm -2 sr" -'A" 1 to be applied to the BFM02 

Jebl values 



atlas. The values derived for the BFM02 target direction 
are Iioo = 0.8 MJy sr _1 and Ay = 0.054 mag, i.e. twice 
as large as the value adopted by BFM02. From these two 
determinations, one obtains an average correction factor 
of 2.9 ± 0.9 by which the BFM02a dust column density 
and DGL intensity estimates have to be multiplied. The 
resulting DGL intensities then become 2.9, 2.3, and 2.3 
xlO- 9 ergs- 1 cm- 2 sr- 1 A- 1 at 3000, 5500, and 8000 A, re- 
spectively. These are to be compared with the Jdgl values 
of 1.0, 0.8, and 0.8 xlO _9 ergs _1 cm _2 sr _1 A _1 as given in 
BFM02a at the three wavelengths. Thus, the DGL inten- 
sity is substantially larger than estimated by BFM02a and 
about equal to their EBL intensities. 

Because the DGL has a similar spectrum as the ZL, 
part of the DGL has been included into the ZL mea- 
surement and therefore has been subtracted from 7 t ot to- 
gether with the ZL. According to BFM02, this part is 
—35 %. Therefore, the remaining DGL correction, to be 
subtracted from Jebl as given in BFM02a, amounts to 
- lxlCrVgs^cm^sr- 1 !- 1 . 

7. THE RESULTING EBL ESTIMATES AND THEIR ERRORS 

My analysis above has revealed several, in part serious 
problems in the BFM02b treatment of the ground-based 
Zodiacal Light observations. 

(1) There is a major puzzle resulting from the compar- 
ison of the BFM02b and BFM98-01 treatment of atmo- 
spheric scattered light: in spite of a large difference in the 
scattered light corrections applied, 7 sca = in BFM98- 
01 vs. 7 sca 24xl(r 9 crgs- 1 cm- 2 sr- 1 A- 1 in BFM02b, 
the derived EBL values are exactly the same to within 
0.1xl0~ 9 ergs^ 1 cm _2 sr _1 A _1 . 

(2) The systematic errors in calibration and atmo- 
spheric scattering, discussed in Sections 4.1-4.3 and 5, 
influence the derived ZL value in the same direction 
(see Table 1); the corrected Izl value at 4650 A is 



smaller than the estimate given in BFM02b by — 14.3- 
19.8xlO _9 ergs~ 1 cm _2 sr _1 A~ 1 . This means that the 
EBL estimates of BFM02 have to be increased by —4.4-6.1, 
-13.7-19.0, and -9.4-13.1 x lO^crgs^cm^sr^A- 1 at 
3000, 5500, and 8000 A , respectively. 

(3) BFM02a have based their estimation of the DGL in- 
tensity on too low a dust column density. The error caused 
by this is — lxlO _9 crgs _1 cm _2 sr _1 A _1 , and is to be sub- 
tracted from the BFM02 Iebl values. This correction is 
small compared with the scattered light and calibration 
corrections described above and does not suffice to com- 
pensate for them (see Table 1); 

(4) The corrected EBL estimates will be in- 
creased to -7.4-9.1, -15.4-20.7, and -10.6-14.3 
xl0- 9 ergs^ 1 cm- 2 sr- 1 A- 1 at 3000, 5500, and 8000 A, 
respectively. These values are 2 to 7 times as high as 
the original BFM02 estimates. Clearly, such high EBL 
values are in conflict with the upper limits of —4.5- 
9xlO _9 ergs -1 cm~ 2 sr~ 1 A~ 1 as derived by Dube, Wickes, 
& Wilkinson (1979), Toller (1983), and Mattila (1990) at 
4000-5100 A. 

(5) BFM02 give for their mean EBL intensities the la 
combined statistical and systematic uncertainties of ±2.5, 
±1.4, and ±1.0xlO- 9 ergs~ 1 cur 2 sr- 1 A- 1 at 3000, 5500, 
and 8000 A , respectively. However, the systematic effects 
discussed in Sect. 3 - 6 of this paper have not been ad- 
equately dealt with in BFM02. The corrections as listed 
in Table 1 are to be considered as systematic errors which 
are not included in the BFM02 error analysis. They can 
be seen to substantially exceed the BFM02 la combined 
statistical and systematic error estimates. 

In summary, I have presented arguments indicating that 
systematic errors of the BFM02 measurement have been 
underestimated. The claim of the paper of a "detection of 
the EBL" appears premature. 



6 



Mattila 



REFERENCES 



Ashburn, E.V. 1954 J. of Atmospheric and Terr. Physics, 5, 83 
Bargcr, A. J., Cowie, L. L., Mushotzky, R. F., & Richards, E. A. 

2001, AJ, 121, 662 
Beckers, J. 1995, in Scientific and Engineering Frontiers for 8-10 m 

telescopes, ed. M. Iye & T. Nishimura (Tokyo: Tokyo Universal 

Academy), 303 

Bernstein, R .A. 1998, The HST/LCO Measurement of the Mean 

Flux of the Extragalactic Background Light (3000-8000A ), 

PhD Thesis, California Institute of Technology (available from 

University Microfilms International) 
Bernstein, R. A. 1999a, in ASP Conf. Ser. 170, The Low Surface 

Brightness Universe, ed. J. I. Davics, C. Impey, & S. Phillipps 

(San Fransisco: ASP), 341 
Bernstein, R. A. 1999b, in ASP Conf. Ser. 193, The Hy-Rcdshift 

Universe: Galaxy Formation and Evolution at High Rcdshift, ed. 

A. J. Bunker & W. J .M. van Bruegel (San Fransisco: ASP), 487 
Bernstein, R. A. 2001, Extragalactic Background Light, in 

Encyclopedia of Astronomy and Astrophysics, ed. P. Murdin 

(Bristol: IOP Publishing and London: Nature Publishing), 778 
Bernstein, R. A., Freedman, W. L., & Madore, B. F. 1999, preprint 
at 4400 A Bernstein, R. A., Freedman, W. L., & Madore, B .F. 2000, 

ApJ, submitted 

Bernstein, R. A., Freedman, W. L., & Madore, B. F. 2002a, ApJ, 

571, 56 (BFM02 Paper I) 
Bernstein, R. A., Freedman, W. L., & Madore, B. F. 2002b, ApJ, 

571, 85 (BFM02 Paper II) 

Bernstein, R. A., Freedman, W. L., & Madore, B. F. 2003, Erratum, 

private communication, communicated via J. P. Huchra, ApJ 

Scientific Editor 
Capaccioli, M., & de Vaucouleurs, G. 1983, ApJS, 52, 465 
Dubc, R. R., Wickes, W. C, & Wilkinson, D. T. 1979, ApJ, 232, 333 
Garstang, R. H. 1991, The Observatory, 111, 239 
Guy, J., Renault, C, Aharonian, F. A., Rivoal, M., & Tavernet, J. -P. 

2000, A&A, 359, 419 
Hartmann, D., & Burton, W. B. 1997, Atlas of Galactic Neutral 

Hydrogen (Cambridge and New York: Cambridge University 

Press) 

Hauscr, M. G. 2001, in IAU Symp. 204, The Extragalactic Infrared 
Background and its Cosmological Implications, ed. M. Harwit & 
M. Hauser (San Fransisco: ASP), 101 

Holtzman, J., et al. 1995, PASP, 107, 156 

Jimenez, R., & Kashlinsky, A. 1999, ApJ, 511, 16 

King, I. R. 1971, PASP, 83, 199 



Krisciunas, K.,& Schaefer, B. 1991 PASP, 103, 1033 

Krist, J., & Burrows, C. 1994 WFPC2 Instrument Science Report 

94-01, The Space telescope Science Institute 
Leinert, C, et al. 1998, A&AS, 127, 1 

Longair, M. 2001, in IAU Symp. 204, The Extragalactic Infrared 
Background and its Cosmological Implications, ed. M. Harwit & 
M. Hauser (San Fransisco: ASP), 505 

Madau, P., & Pozzetti, L. 2000, MNRAS, 312, L9 

Mackie, G. 1992, ApJ, 400, 65 

Mattila, K. 1990, in IAU Symp. 139, The Galactic and Extragalactic 
Background Radiation, ed. S. Bowyer & Ch. Leinert (Dordrecht: 
Kluwer), 257 

McClatcey, R. A., Fenn, R. W., Sclby, J. E. A, Volz, F. E., & Garing, 
J. S., 1978, in Handbook of Optics, ed. W. G. Driscoll & W. 
Vaughan (New York: McGraw-Hill), 14-14 

Middlcmass, D., Clegg, R. E. S., & Walsh, J. R. 1989, MNRAS, 239, 
5P 

Pagcl, B. E. J. 2002, in Chemical Enrichment of Intracluster and 
Intergalactic Medium, in ASP Conf. Ser. 253 , ed. R. Fusco- 
Femiano & F. Matteucci (San Fransisco: ASP), 489 

Peebles, P. J. E. 2001, in IAU Symp. 204, The Extragalactic Infrared 
Background and its Cosmological Implications, ed. M. Harwit & 
M. Hauser (San Fransisco: ASP), 5 

Piccirillo, J. 1973, PASP, 85, 278 

Pozzetti, L., & Madau, P. 2001, in IAU Symp. 204, The Extragalactic 
Infrared Background and its Cosmological Implications, ed. M. 
Harwit & M. Hauser (San Fransisco: ASP), 71 

Primack, J. R., Bullock, J. S., Somerville, R. S., & MacMinn, D. 
1999, Astroparticle Physics, 11, 93 

Racine, R. 1996, PASP, 108, 699 

Renault, C, Barrau, A., Lagache, G., & Pugct, J.-L. 2001, A&A, 
371, 771 

Roach, F.E., & Megill 1961, ApJ, 133, 228 
Roddier, F. 1995, Ap&SS, 223, 109 

Schlcgcl, D. J., Finkbeincr, D. P., & Davis, M. 1998, ApJ, 500, 525 
Surma, P., Seifert, W., & Bender, R. 1990, A&A, 238, 67 
Staude, H. J. 1975, A&A, 39, 325 

Toller, G. N. 1981, Ph.D. Thesis, State University of New York 
Toller, G. N. 1983, ApJ, 266, L79 

Uson, J. M., Boughn, S. P., & Kuhn, J. R. 1991, ApJ, 369, 46 
Wheelock, S. I., Gautier, T. N., Chillcmi, J., et al. 1994, IRAS 

Sky Survey Atlas Explanatory Supplement, JPL Publication 94-11 

(Pasadena: JPL) 
Wright, E. L. 2001, ApJ, 553, 538 



